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Abstract 

The effects of initial b quark bound state for the semi-inclusive decays B —* 
K(K*)X are studied using light cone expansion and heavy quark effective 

a.: 

theory methods. We find that the initial bound state effects on the branching 



ratios and CP asymmetries are small. In the light cone expansion approach, 
the CP-averaged branching ratios are increased by about 2% with respect to 
the free 6-quark decay. For B° -> K~(K*~)X, the CP -averaged branching 
ratios are sensitive to the phase 7 and the CP asymmetry can be as large as 
7% (14%), whereas for B~ — > K°(K*°)X the CP-averaged branching ratios 
are not sensitive to 7 and the CP asymmetries are small (< 1%). The CP- 
averaged branching ratios are predicted to be in the ranges (0.53 ~ 1.5) x 10 -4 
[(0.25 ~ 2.0) x 10~ 4 ] for B° -> R-{K*-)X and (0.77 ~ 0.84) x 10" 4 [(0.67 ~ 
0.74) x 10~ 4 ] for B- -> K °(K*°)X, depending on the value of the CP violating 
phase 7. In the heavy quark effective theory approach, we find that the 
branching ratios are decreased by about 10% and the CP asymmetries are 
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not affected. These predictions can be tested in the near future. 
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I. INTRODUCTION 



There have been considerable experimental and theoretical efforts to understand the 
properties of B decays. These studies have provided important information about the mech- 
anism for B decays and the origin of CP violation. In the next few years large quantities 
of experimental data on B decays will become available. It is hoped that one will obtain 
even more important information in understanding the mechanism for B decays and the 
mechanism for CP violation. In particular, charmless hadronic B decays have played an 
important role in the determination of the CP violating parameter 7 in the Standard Model 
(SM) fp-fl. While most of the studies have concentrated on the exclusive B decay modes 
for CP violation, there are also some studies for semi-inclusive decays |2|||. At the quark 
level the relevant Hamiltonian for B decays in the SM is well understood. The major un- 
certainties for these decays come from our insufficient understanding of the long distance 
strong interaction dynamics involved in these decays. There are several methods which have 
been used to estimate the decay amplitudes, including naive factorization, QCD improved 
factorization and methods based on symmetry considerations. 

Recently it has been argued that in the heavy quark limit, factorization is a good approx- 
imation Q and several processes have been calculated ||. Leading QCD corrections to the 
naive factorization can be studied for exclusive decays in a systematic way. In the calculation 
of exclusive decays, the hadronic matrix elements can be factorized and strong interaction 
dynamics can be parameterized into the relevant decay constants, light cone distribution 
amplitudes and transition form factors. At the present time, the light cone distribution am- 
plitudes and transition form factors are not well known which introduce uncertainties in the 
calculations. Of course one should keep in mind that there may be large corrections of order 
^QCDj^b which needs further study. From quark hadron duality consideration, inclusive 
decays can be represented by quark level calculations and the uncertainties may be small. 
It is believed that theoretical calculations for exclusive decays contain more uncertainties 
than inclusive decays. Of course when going completely inclusive, there are less information 
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that can be extracted about strong and weak interaction dynamics and CP violation, and 
it is experimentally hard to identify final states inclusively. In this paper we will take the 
way in between by studying semi-inclusive decays following Ref. [0] in the hope that one 
may be able to reduce some of the hadronic uncertainties in exclusive decays on one hand, 
and still be able to obtain important information about B decays and CP violation with 
clear experimental signal on the other. We will study the charmless semi-inclusive decays 
B — > KX and B — > K*X. Here the X indicates states containing no charmed particles. 

The decay modes B — > K(K*)X have been studied before [@|§. 111 previous studies, 
several effects were treated phenomenologically, such as the number of colors was taken as 
an effective number and treated as a free parameter, the gluon virtuality q 2 in the penguin 
diagrams was assumed to be around ml/ 2, and the bound state effects of b-quark inside 
the B meson was modeled by assuming its momentum to obey a Gaussian distribution. 
To have a better understanding of these decays, it is necessary to carry out calculations in 
such a way that the phenomenological treatments can be improved with better theoretical 
understanding. It has recently been shown that it is indeed possible in the heavy quark 
limit to handle most of the problems in exclusive B to two light meson decays from QCD 
calculations We will use the same formalism in our study of semi-inclusive decays in the 
factorization approximation, paying particular attention to the initial bound state effects. 

The problems treated in the case of exclusive decays are different in some ways from the 
semi-inclusive decays studied here. The problems associated with the number of colors and 
the gluon virtuality can be treated the same way, but the initial b quark bound effects in 
semi-inclusive decays arise in different form from those in exclusive decays. In the exclusive 
decay case, the b quark bound state effects are taken care by decay constants and transition 
form factors. In the semi-inclusive case, there are contributions which, in the free quark 
decay approximation, can be viewed as a b quark decay into a meson and another quark. 
One needs to treat initial b quark bound state effects on more theoretical ground. This will 
be the main focus of this paper. We will study this problem using two different methods 
with one based on light cone expansion and another based on heavy quark effective theory. 
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To further reduce possible uncertainties associated with form factors, we will choose 
processes which have the least numbers of hadronic parameters beside the ones related to 
the initial bound state effects. We find that the following processes are particularly good 
for this purpose, 

B° -> K~X, B~ -> K°X, 

B° -> K*-X, B~ -> K*°X. (1) 

For these processes, the transition form factors for B — > K and i? — > K* do not show 
up in the factorization approximation because the bi-quark operator sTb (here T is some 
appropriate Dirac matrices) does not change the electric charge of the initial particle B and 
the final particle K(K*). Therefore for these processes there are only the K(K*) decay 
constants and parameters related to the initial bound state effects if small annihilation 
contributions are neglected. 

The paper is arranged as follows. In Section II, we will study the decay amplitudes in 
the SM for the semi-inclusive B — > K(K*)X decays. In Section III, we will study the light 
cone and heavy quark effective theory formulation of the initial bound state effects on these 
semi-inclusive decays. And in Section IV, we will carry out numerical analyses of the energy 
spectra of the K(K*), branching ratios and CP asymmetries in B — > K(K*)X, and draw 
our conclusions. 



II. DECAY AMPLITUDES IN THE HEAVY QUARK LIMIT 

In this section we study the short distance decay amplitudes for semi-inclusive B — > 
K(K*)X decays. The effective Hamiltonian for charmless B decays with AS = 1 at the 
quark level is given by 

H eff = ^ {^^(ciOi + c 2 2 + CnO n ) + V cb V c * s ^ C n O n \ . (2) 

I n=3 n=3 J 

Here O n are quark and gluon operators and are given by 
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Ox = (SiUj)v-A{Ujbi) V -A, 2 = (SiUi) V -A{Ujbj)v-A, 

0-3(5) = (Sibi)v-A^2((fjqj)v-{+)A, 04(6) = (Sibj)v-A XX^j'^O V-(+)A, 
q> q> 

3 3 
°7(9) = 2{si h i)v-AY. e <i'% < l'j)v+{~)A, Cs(io) = 2(^ & i)v-A^ e 9'(^%)v+(-)A, 

On = ^ro 6 *^G^(l + 75)&„ (3) 

where (V ± A)(V ± A) = 7 M (1 ± 75)7^(1 ± 75), q' = u,d,s,c,b, e q i is the electric charge 
number of the q' quark, A a is the color SU(3) Gell-Mann matrix, i and j are color indices, 
and G^ u is the gluon field strength. 

The Wilson coefficients c n have been calculated in different schemes || . In this paper we 
will use consistently the NDR scheme. The values of c n at \x ~ m& with the next-to-leading 
order (NLO) QCD corrections are given by || 

Cl = -0.185, c 2 = 1.082, c 3 = 0.014, c 4 = -0.035, c 5 = 0.009, c 6 = -0.041, 

c 7 = -0.002a em , c 8 = 0.054a em , c 9 = -1.292a em , c w = -0.263a em , c u = -0.143. 

Here a em = 1/137 is the electromagnetic fine structure constant. 

In order to make sure that the observed events are from rare charmless B decays, and 
other processes, such as B — > D(D*)X f — > K(K*)X" , do not contaminate the direct rare 
decay of B — > K(K*)X due to short distance interaction, we will make a cut on the K(K*) 
energy which will be set at Ek,k* > 2.1 GeV. It has been shown that this cut can eliminate 
most of the unwanted events while leave most of the events induced by short distance 
contributions @ because the matrix elements of the type < K(K*)\ji\0 >< X\j 2 \B > 
would results in a fast K(K*) in the final state. The resulting events will resemble two 
body type of decays with one of them be the K{K*) and another, back-to-back against the 
K(K*), will be X with small invariant mass M\. With the cut E K>K * > 2.1 GeV, M\ < 5.7 
GeV 2 . 

The hadronic matrix element for a specific operator < XK\0\B > is difficult to calculate 
at present. We will use factorization approximation to estimate it. The factorization ap- 
proximation has been shown to hold in the heavy quark limit for exclusive B decays into two 
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light hadrons. The leading contribution for an operator which can be written as a product 
of two currents ji = sTiq' and j 2 = q'T 2 b with I\ carrying appropriate Lorentz and Dirac 
indices, O = j \ ■ j 2 , is given by 

< XK\0\B > fact = < K\jt\0 >< X\j 2 \B > + < X\j[\0 >< K\f 2 \B > 

+ <XK\j 1 \0><0\j 2 \B>. (4) 

The second term on the right-hand-side in the above represents the Fierz transformed fac- 
torization terms with j[ = q'T[q' and j 2 = sT 2 b. The third term is usually referred to as the 
annihilation contribution. 

B — > KX is a many-body decay, which is different from two-body decays. There are 
more ways of factorization for a many-body decay, such as < Xii^|ji|0 >< X^^l-B > and 

< X 2 |j(|0 >< X' 2 K\j' 2 \B >, with X = X x + X[ = X 2 + X' 2 . The three terms in Eq. (§ 
corresponding to the cases: < X±\ =< 0|, < X 2 \ =< 0| and < X[\ =< 0|, respectively. 
For B — > KX with a cut Ex > 2.1 GeV, the final state X has a small invariant mass. 
This is a quasi-two-body decay, with K and X moving rapidly apart in opposite directions. 
The probability of forming the final state < X\K\ with < X\\ ^< 0| is less than the 
probability of forming the simple final state < K\. This suggests that the contribution of 
the configuration < XiKlj^O >< X[\j 2 \B > is dominated by < i^|ji|0 >< X\j 2 \B >. 
Likewise, the contribution of the configuration < X 2 |jj|0 >< X 2 K\f 2 \B > is dominated by 

< X|jj|0 >< K\j' 2 \B >. The cases with \Xi > and \X 2 > not equal to |0 > are also higher 
order in a s and therefore a s power suppressed. We will neglect them in our later discussions 
which also eliminate the third term in Eq. (|]). 

The above approximation is also supported by explicit calculation of the bremsstrahlung 
process, b — > Kq'g, which represents some of the a s order corrections. It has been shown, 
in a similar situation of b — > <ps and b — > <psg, that the bremsstrahlung contributes less than 
3% of the total branching ratio 0] . One can easily obtain from Ref . |7j an estimate of the 
contribution for the processes considered here. The bremsstrahlung contribution is small. 
Eq. will adequately approximate the leading contributions and we will work with this 
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approximation. 

In the heavy quark limit, a class of radiative corrections in powers of a s , which does 
not change the form of the operators, can be included for the matrix elements. For a local 
operator the correction can be parameterized as the following, similar to the exclusive decays 
discussed in Refs. 

oo 

< XK\0\B >=< XK\0\B > fact [l + Y. r n< + 0(A QCD /m b )], (5) 

n=l 

where < XK \0\B > f act denotes the naive factorization result. Aqcd ~ 0.3 GeV is the strong 
interaction scale. The second and third terms in the square bracket indicate, respectively, 
higher order a s and Agc^/m;, corrections to the factorized matrix element. 

Similar arguments can be made for B —>■ K*X decays also. For the < K(K*)\ji\0 >< 
X\j 2 \B > type, the decay amplitudes involves the K(K*) decay constants, while for the 

< X|jj|0 >< K(K*)\j' 2 \B > type, it involves the transition form factors from B to K(K*), 
and the < XK(K*)\ji\0 >< 0|j2|-£> > type involves the B decay constant. 

If all three terms in Eq. (|j) contribute with the same order of magnitude, the accumulated 
uncertainties will be substantial due to large uncertainties in the transition form factors and 
the B decay constant. Fortunately we find that for B° -> K~{K*~)X and B~ -> K°(K*°)X, 
only the first and the third types of terms in Eq. contribute due to electric charge 
conservation. This eliminates possible uncertainties from the transition form factors. Also 
as argued before the third term can be neglected because it is subleading and a s power 
suppressed. There is only one term present, which considerably simplifies the calculation. 

Using the effective Hamiltonian in Eq. (0), we obtain 

A(B - KX) =& £ V qb V*J K [A«P£ < X\q%(l - l5 )b\B > 

V * q=u,c 

+ B"<X\q\l- l5 )b\B>l 
A(B - K*X) = ^ E V qh V* 8 m K *f K *A^ < X|^(l - l5 )b\B >, (6) 

V ^ q=u,c 

where q' = u and d for B° and B~, respectively. The decay constants are defined as 

< K\s^(l -75)g'|0 >= if K PR and < K*(X)\s^(l - j 5 )q'\0 >= m K *f K *e^. We adopt the 
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standard covariant normalization < B\B >= 2Eb{2it) 3 5 3 (0). The coefficients A q (A q ) and 
B q are given by, for B° -> K~(K*-)X 



A q (A q ) = a\ + a\ + a q 10 + a 



10a > 



B q = (ai +a q R + ai 



2m 



K~ 



m u + m s 



For B- -» K°(K*°)X } 



A q (A q ) — a\ — -a? + a? 0a , 



2 m d + m s 

Including the lowest a s order corrections in Eq. (^j), af are given by 

u , ci ajCj 

a$ = 0, 

a! = c 4 + ^ + ^% [c 3 (F P + G P (s s ) + G P (s 6 )) + Cl G P (s q ) 

iV 47T iV 
+ (C 4 + C 6 ) G p( s f) + CnGp.ll 

«e = c 6 + + [c 3 (G P (s s ) + G P (s 6 )) + ciG' P (s q 

b 

+ (c 4 + c 6 ) ^ G' P (sf) + ciiG P11 
/=« 

9 i C 7 

8 = ° N' 



3 b 3 
(c 8 + cio)- 51 e f G' P (s f ) + c 9 -(e s G P (s s . 



ebGp(sfc)) 



a Cg tt s Gp 



q a s Cp 
ai0a - 4^ AT 



3 3 
(c 8 + ci )~ e f Gp(s f ) + c 9 -(e s G P (s s ) + e b G P (s b )) 

f=u 



(7) 



(8) 



(9) 



where N = 3 is the number of colors, Cp = (N 2 — 1)/(2N), and Sf = m 2 /m 2 . The other 
items are given by 



F P = -12 In - 18 + fL 
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Jo 1 — X 




G' K ,(s) = 0, 



G 



K* ,11 ~ 



(10) 



where the subscript P can be K or K*, indicating that the coefficients af are process de- 
pendent. <Pk{x) and <f)°x{x) are the twist-2 and twist-3 kaon meson distribution amplitudes, 
respectively, (fix* (%) is the leading twist distribution amplitude for the longitudinally polar- 
ized K*. In this paper we will take the following forms for them ||, 



limit. The number of colors should not be treated as an effective number, but has to be 
3 from QCD. The results are, in principle, renormalization scale and scheme independent. 
The problem associated with the gluon virtuality k 2 = (1 — x)m 2 B in the naive factorization 
calculation is also meaningfully treated by convoluting the x-dependence with the meson 
distribution amplitudes in the functions Gp(s) and G' P (s). 



In this section we study the decay rates for B — > K(K*)X, taking into account b quark 
bound state effects, using two different methods, the light cone expansion method and the 
heavy quark effective theory method. 

We will work out, in detail, the formulation for B — > KX in the following. The results 
for B — > K*X can be easily obtained in a similar way. Without taking into account the 



<pK,K*{.x) = 6x(l - X), (f> K (x) = 1. 



(11) 



The amplitudes in Eq. (|6|) are from perturbative QCD calculation in the heavy quark 



III. INITIAL BOUND STATE EFFECTS 
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initial bound state effects, that is in the free b quark decay approximation, the decay can 
be viewed as the two body process b — > Kq' and one obtains [Q] 

T(B -> KX) « r(6 -> Kg') = ^K|a| 2 + |/3| 2 )m 6 , 

« = % E = % E *W* tf - (12) 

V ^ q=u,c V ^ q=u,c 

If the b quark mass is infinitively large, Br(B — > is equal to Pr(6 — > K(K*)q'). 

However due to initial b quark bound state effects there are corrections . We now proceed 
to study the initial bound state effects on the decay rates. 

The differential decay rate for B — > KX in the B rest frame, following the procedure in 
Ref. H, is given by 

i r/ 3 P 

^ KX) = ^ bW ^- k E^Wft " ^ " Px)|A(fl - KX)| 2 . (13) 
Using / d 4 y exp[-zy ■ (P B - P K - P x )} = {2tx) a 5\P b - P K - P x ), we have 
£(2tt)V(P b - P K - P X )\A{B - KX)| 2 = f K J2 I d "y e-^- 13 *-^ 

X X J 

x [\a\ 2 P^ K < B\h v {l - ls)q'\X >< X\q\(l - l5 )b\B > 
+ \P\ 2 < B\b(l + l5 )q'\X >< X\q'(l - l5 )b\B >] 

= fl\ A e^ P H\a\ 2 P^ < B\[jt(0),Uy)]\B > +|/3| 2 < B|[jt(0), J(y)]\B >), (14) 

where j M = q'^f^l - 7s)& and J = q'(l - 75)6. 

Computing the current commutators one obtains 

£(27r) 4 a 4 (p B - - Px)|a(p - kx)| 2 

X 

J d 4 y e*y p «id a A q ,(y)] < B\b(0)^(l- l5 )U(0,y)b(y)\B > . (15) 



x 



In the above we have assumed m q > = and used 

W(x), q'(y)} = 1(7 ■ d)iA q ,(x - y)U(x, y), (16) 

with 
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U(x,y) = Vexp[ig s f dz^G^z)], 
J y 

Mv) = ~j^y 3 J d 4 k e-^e(k°)6(k 2 ), (17) 

where U(x, y) is the Wilson link, G M is the background gluon field, and e(x) satisfies e(|x|) = 1 
and e(— \x\) = —1. 

The matrix element < i?|6(0)7 /3 (l — j5)U(0, y)b(y)\B > which is equal to < 
B\b(0)^U(0,y)b(y)\B > from parity consideration contains all information about initial 
bound state corrections. It is, however, difficult to completely evaluate it due to non- 
perturbative effects. In the following we attempt two calculations: one using light cone 
expansion, and the other using heavy quark effective theory. 

A. Light Cone Expansion Estimates 

In general one can decompose the matrix element, < B \ b(0)j^U(0, y)b(y)\B >, in the 
following form 

< B|6(0)/E7(0, y)b(y)\B >= 2[P^F(y 2 , y ■ P B ) + y^G(y 2 , y ■ P B % (18) 

where F(y 2 ,y ■ Pb) and G(y 2 ,y ■ Pb) are functions of the two independent Lorentz scalars, 
y 2 and y ■ P B - 

Since we are interested in having large kaon energy > 2.1 GeV and small invari- 
ant mass for the X, the dominant contribution to the y integration in Eq. ( |T5D will be 
from the light cone region y 2 < \jE\, which suggestes that, as a good approximation, 
< i?|6(0)7 /3 ?7(0, y)b{y) \B >^ 2P B F(0,y ■ Pb)- This approximation is also supported by 
the fact that the function A q /(y) has a singularity at y 2 = while away from light cone it 
vanishes. Carrying out a Fourier transformation ||, 

F(0,yP B ) = jdO ''"•'••'/(C). (19) 

and inserting the above into Eq. (|15D, we arrive at 
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£(27t) 4 5 4 (Pb -Pk- Px)\A{B - KX) 
x 

= Mf K {2\a\ 2 P%P K -P B + \P\ 2 P%) 



x 



J di 5[{iP B - P K ) 2 MP Ba - P Ka )f(0- (20) 



We finally obtain the decay distribution as a function of Ek 

dT(B^KX) f 2 K . , , 2 . , |<a|2 v p ,.2£ 



(4|a| 2 P| + |/3| 2 )^ / (^). (21) 

\ ' 1 11 r> 



dEx 2nmB ^ ' tub 

Carrying out similar calculations, we obtain the differential decay rate for the B — > K*X 

decay 

K " X) - ^4|«/(^), (22) 
dEx* ZirmB mg 

where a* = (G F /V2) J2 q =u, c V qb V* s Ai . 

It is interesting to note that the same distribution function /(£) appears in both B — > KX 
and P> — ► i^*X cases. It is also interesting to note that, in the approximation made in this 
section, the function /(£) is the same as that in B — > X7 || and semi-leptonic decays 
B — > X/i/ |K|. These decays have been studied in details. Experiments in the future 
will measure the differential distributions for these decays and, therefore, provide detailed 
information about /(£). We can use this information in the calculation to reduce error. 
One may also turn the argument around to use the decay modes discussed here to provide 
constraints on the form of the distribution function /(£). Before the detailed experimental 
information becomes available, we have to make some theoretical modeling for our numerical 
analysis which will be discussed later. 



B. Heavy Quark Effective Theory Estimates 

We note that the simple expressions for the decay distributions in Eqs. ([H]) and flZ2"|) 
hold to leading order in light cone expansion. When higher order contributions are included 
the expressions will not be so simple. To have some idea about the sensitivity of the results 
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to other corrections, in the following we also estimate the corrections to the free b quark 
decay rates using heavy quark effective theory. 

If the b quark is heavy, the decay products all have large energy and to a good approx- 
imation can be treated as free quarks. In that case, {7(0, y) ~ 1 because the background 
gluon field can be approximated to vanish, we then have 

< B\b(0)^U(0,y)b(y)\B >^< B\b(0)^b(y)\B > . (23) 

If the b quark is infinitively heavy, the above matrix element is simply given by 2P^e~ tmbV y , 
where v is the four velocity of the B meson satisfying v 2 — 1. Since the b quark has finite 
mass, there will be corrections. We now estimate the leading 1 / m\ corrections following the 
procedure outlined in Ref. ||. In the heavy quark effective theory, the b(x) quark field can 
be expanded as 

b(x) = e- im " v - x {l + i 1 ■ D T /(2m b ) + v ■ D 1 ■ D T /(4m 2 b ) - ( 7 • D T f /(8m 2 b )}h(x) + 0(l/m 3 b ) 

+ (terms for anti-quark), 
D% = D" - v^v ■ D, 

= d» -ig s G»{x). (24) 
Using the above expressions and keeping l/m& terms, we obtain [§] 

< B\b(0)^b(y)\B > = 2m B e- im «{^ - ^-(V + « ■ V^){£ ' 

-\{y 2 -{vy)VA)}, (25) 

where 

1 



A = T~ < B\hg s G^h\B >, 

l4 = -t^— < B\h(tD T ) 2 h\B > . (26) 
2m b 



We note that the expansion in Eq. fl25|) is different from the light cone expansion as can be 
seen from the above expression that some y 2 terms are kept. The expansion is truncated at 
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order l/m& in Eq. (p5|). The truncation of the 1/mb expansion enforces the use of the quark 
level phase space, instead of the hadron level phase space. 
Inserting the above expression into Eq. ([13]), we have 

F(B - KX) « f m 6 [l««(l + ^4 " f 4) + l/?| 2 (l - ^ + ( 2 ?) 

In the approximation made here, the distribution of E% is a delta function with the peak 
at E K = rrib/2. 

Carrying out similar calculations, we obtain the decay rate for the B — > K*X decay, 

,2 



r(fi - « ^m 6 K| 2 m 2 (l + 7 ^ - f ^|). (28) 

87r 6mj 6 



It is clear that in the limit of large m&, that is /i 2 5 /m 2 — > 0, the result reduces to the 
free b quark decay b — > K(K*)q' result as expected. 

The expressions for the decay rates, in the approximation we are working with, are 
simple, allowing easy analysis. In the case of the light cone expansion method, one needs to 
have detailed knowdege of distribution function /(£) for numerical analysis. Although the 
detailed shape is not known, we do know some properties [PUlQfl. When integrating £ from 
to 1, Jq 1 must give 1 due to current conservation. If the decay can be considered to 

be a free b quark decay, then U(0,y) = 1 because no background gluon field exists, and the 
b quark field is given by b(y) = e~ iy ' Pb b(0), one obtains 

/(£)=*(*-—)■ (29) 
m B 

We can also estimate the mean < £ >= Jq 1 and the variance a 2 = Jq 1 g?££ 2 /(£) — 

< £ > 2 using heavy quark effective theory. They are given by [|l0|,p| 



<t>= — [i+Aa4-/a 

m B bmf 



The small value for a 2 implies that the distribution function is sharply peaked around 
m b /m B . 

15 



To go further we take the following parameterization for the distribution function (10 



f(o= v-y+w <3i) 

where N is a normalization constant which guarantees Jq = 1. This function reduces 

to a 5-function with the peak at a as 6 — > 0. Comparing with Eq. (|29"D, in this limit 
a = rrib/mB- Once the parameters c and <i are given, the parameters a and 6 can be fixed 
by comparing with < £ > and cr 2 . Unfortunately we do not know the values for c and d at 
present. We will take c and d to be free parameters and vary them to see how the energy 
spectra of K(K*), branching ratios and CP asymmetries are changed. 

IV. RESULTS AND DISCUSSIONS 

We are now ready to present our numerical analysis. We will make theoretical predic- 
tions for the kaon energy spectra dT /dE K (*), CP-averaged branching ratios and direct CP 
asymmetries defined as 

Br ave (B -> K(K*)X) = ^[Br(B -> K(K*)X) + Br(B -> K(K*)X)\, 
A (R^K(K*\Y\ r(g -> K(K*)X) - T(B -> g(g^g) 

A CP (i? - /f (if )x) = r{B ^ K(K * )x) + r{ B^K(K*)xy (32) 

For the numerical analysis, we need to know the values for the parameters involved. 
Some of them are well determined. In our numerical calculations we will use the following 



values for the relevant parameters ||11|| : mj = 4.9 GeV, m c = 1.5 GeV, m s = 120 MeV, 
m d = 4 MeV, m u = 2 MeV, \V US \ = 0.2196, \V cb \ = 0.0402, \V ub /V cb \ = 0.085, f K = 160 
MeV, f K * = 214 MeV, a s (M z ) = 0.118. We keep the CP violating phase 7 to be a free 
parameter and vary it to see how the branching ratios and CP asymmetries depend on it. 
The HQET parameter /z 2 , can be extracted from the B* — B mass splitting: /i 2 = 3(m 2 3 » — 



in 



B 



)/4 ~ 0.36 GeV 2 , while /1 2 is less determined. A calculation of QCD sum rules gives 



/1 2 = (0.5 ± 0.2) GeV 2 @, which is consistent with /x 2 = (0.45 ± 0.12) GeV 2 from a recent 



lattice QCD calculation [D^]. We will use /x 2 = 0.5 GeV 2 for our numerical calculations. 
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In the case of light cone expansion, we also need to specify the distribution function 
/(£). We will assume it to be the form given in Eq. (plf). To have some idea how the 
kaon energy spectra, branching ratios, and CP asymmetries depend on the form of the 
distribution function, we consider two very different forms [|14]]: (i) preset c = d — 1, in that 
case a = 0.9548 and b = 0.005444 determined by the known mean value and variance of 
the distribution function; (ii) preset c = d = 2, in that case a = 0.9864 and b = 0.02557 
determined by the same mean value and variance of the distribution function. 

In Figs. 1-4, we show the kaon energy spectra in B — > K(K*)X decays computed in the 
light cone expansion approach, assuming 7 = 60°. The solid and dashed curves correspond, 
respectively, to the parameter set (i) and (ii) for the distribution function. The kaon energy 
spectra are a discrete line at E K (*) = m&/ 2 in free b quark decay approximation, which is 
not shown in the figures. We see that initial bound state effects stretch the spectra over the 
full kinematic range < E K {*) < rriB/2 and the kaon energy spectra depend strongly on the 
form of the distribution function. However, we note that all the spectra have more than 97% 
of events with E K <*) > 2.1 GeV. This implies that if the integrated branching ratios and CP 
asymmetries are measured with E K (*) > 2.1 GeV, the effects from the detailed shape of the 
distribution function are small. 

We show the CP- averaged branching ratios, in Figs. 5-8, and the CP asymmetries, in 
Figs. 9-12, in B — ► K(K*)X as a function of the CP violating phase 7. The solid curves 
are the results from the light cone expansion using the parameter set (i) for the distribution 
function, while the dashed curves are from the free b quark decay approximation. The initial 
bound state effects encoded in the distribution function almost cancel completely in the CP 
asymmetries in B — > K*X, so that the solid and dashed curves coincide in Figs. 11 and 
12. We find that the shifts in the branching ratios and CP asymmetries are negligible if the 
parameter set (ii) instead of (i) for the distribution function is used, indicating that both 
the branching ratios and the CP asymmetries are insensitive to the detailed shape of the 
distribution function. 

One can clearly see from Figs. 5-8 that the differences between the solid and dashed 
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curves are small, about 2%. This implies that according to light cone expansion estimates, 
the initial bound state effects increase the CP-averaged branching ratios for B — > K(K*)X 
by about 2%, largely because the B —>■ K(K*)X phase space is used, which is larger than 
the b — > K(K*)q' phase space used in the free 6-quark and heavy quark effective theory 
calculations. The branching ratios for B° — > K~(K*~)X are sensitive to 7, varying from 
0.53(0.25) x 10~ 4 to 1.5(2.0) x 10" 4 , whereas the branching ratios for B~ -> i?°(K*°)X 
are not sensitive to 7, varying from 0.77(0.67) x 10~ 4 to 0.84(0.74) x 10~ 4 . The above 
sensitivities to 7 can be easily understood by noticing that the tree operators contribute 
to B° — > K~(K*~)X decays but not to B~ — > K°(K*°)X decays when small annihilation 
contributions are neglected, resulting in strong dependence on V u bV* s for the former, but 
not for the latter. 

For the same reasons, the CP asymmetries are expected to be much larger in B° — > 
K~ [K*~)X than in B~ — > K°(K*°)X . The differences between the solid curves and dashed 
curves in Figs. 9 and 10 are very small, about 1%. This implies that according to light 
cone expansion estimates, the initial bound state effects increase the CP asymmetries in 
B — > KX by about 1%. They do not affect the CP asymmetries in B — > K*X. The CP 
asymmetries in B° — > K~ (K*~)X can be as large as 7%(14%), but very small (< 1%) in 
B~ -> K°(K*°)X, as expected. 

The heavy quark effective theory estimates of the initial bound state effects are always 
to reduce the branching ratios at the level of 10% as can be seen from Eqs. ([27]) and (|28|) if 
(/% = n 2 g is used. In fact within the allowed range for //^ the initial bound state effects tend 
to reduce the branching ratios. The CP asymmetries are the same as those obtained by free 
b quark decay approximation. 

The three estimates (free quark decay approximation, light cone expansion and heavy 
quark effective theory method) carried out here all give the same order of magnitudes for 
the branching ratios and CP asymmetries which are also the same order of magnitudes as 
those obtained in Ref. 0. The initial bound state effects are at the order of 10% of the 
free b quark decay estimates. The differences between different methods may be viewed as 
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uncertainties in the estimates. The branching ratios are of order 1CT 4 and are within the 
reach of the B factories. The CP asymmetries in the neutral B modes -B — > K~(K*~)X 
are large and can be measured at the B factories. When more data become available, one 
may obtain interesting information about hadronic effects and also information about the 
CP violating phase 7. 
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FIG. 1. Kaon energy spectrum in B° — > A" X. In Figs. 1-4, the solid curves are for (i) 
c = d = 1; the dashed curves are for (ii) c = <i = 2. 
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FIG. 5. CP-averaged branching ratio for B° — > In Figs. 5-10, the solid curves are 

for the ligh cone expansion with (i) c = d = 1; the dashed curves are for the free 6 quark decay 
approximation. 
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FIG. 6. CP-averaged branching ratio for B — ► 
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FIG. 9. CP asymmetry in B° -► K~X. 
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FIG. 10. CP asymmetry in B~ -»• K°X. 
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